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Abstract 

We discuss the effects of QCD corrections to the on-shell decay t — > bW. We resolve 
the scale ambiguity using the Brodsky-Lepage-Mackenzie scheme, and find that the appro- 
priate coupling constant is a^ s (0.122mt). The largest long distance contribution comes 
from the definition of the on-shell mass of the top quark. We note that QCD corrections 
to the electroweak p parameter are extremely small when the p parameter is expressed in 
terms of the top quark width. 



There has been much discussion in the literature of QCD corrections to electroweak 
processes involving the top quarkQ' nS i. One of the motivations for this has been an 
attempt to extract as much information as possible from electroweak radiative corrections 
to measurements such as T(Z — > bb) and the p parameter. With the measurement of 
possible signals of top production at the Tevatron^ and the observation of electroweak 
radiative corrections at LEP0, there is now the exciting possibility of using a measured 
value of m t to extract information about the Higgs boson or new physics beyond the 
Standard Mode from precise electroweak data. 

In a previous paper @, we have discussed QCD effects on electroweak corrections, 
specifically the p parameter. We showed that the QCD corrections were dominated by 
perturbative effects involving momentum transfers on the order of m t . Since any top 
quarks appearing in radiative loops have a high virtuality, there are no sizable tt thresh- 
old corrections. The largest nonperturbative effect is associated with expressing the p 
parameter through the top quark mass measured at the pole, m t . This is an uncertainty 
that arises when a quantity which is well defined at short distances is written in terms 
of a long distance parameters and brings about corrections on the order of KqcD/ m t Q- 
When the p parameter is written in terms of quantities defined at short distances, all 
nonperturbative effects are suppressed by a factor of {^QCD/mt) 4 - 

In this paper, we resolve all normalization point ambiguities in the top decay rate in 
the order a s by using the method of Brodsky, Lepage, and MackenzieS. The basic idea 
of the BLM prescription is that when one changes the normalization point from a s (pi) 
to a s {p2)-, the amplitude will gain a higher order term proportional to the number of 
light quark flavors. When a result is expressed in two different scales, the difference is 
an rif dependent higher order in a s term. A poor choice of scale will contain a large rif 
dependent higher order corrections. If one wants to find an appropriate scale to use at 
order a s , one should calculate the rif dependent part of the a 2 correction. There exists 
a scale for which the rif dependent a 2 correction is zero. This scale roughly corresponds 
to the weighted average of gluon momenta. It was argued by BLM that this scale is the 
physically relevant scale for any given problem. 

In analyzing QCD corrections to the p parameter, we found that the relevant cou- 
pling constant according to the BLM prescription was af IS (0.154 : mt). We will proceed to 
analyze the QCD corrections for the on- shell decay rate of t — > bW, which is by far the 

1 The intrinsic ambiguity of order Kqcd in the definition of the quark on-shell mass was also recently 
emphasized in detail by Bigi et. al. . 
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dominant decay process for the top in the Standard Model. All calculations are performed 
in the limit m? >> mL. The leading in a H correction can be found in the literature for 
finite ,J .In al, calculations, the Cabibbo-Kobyasfii-Maskawa matrix element V a is taken 
to be unity. In this limit, the width is given to the leading order in a s by! 



where T = ^f^ - is the tree-level width. We wish to resolve the scale ambiguity in this 
correction by calculating the n/ dependent part of the 0(a 2 ) correction, and quantify this 
scale by using the BLM prescription. 

The rif dependent virtual correction can be found by considering all one gluon contri- 
butions with an additional vacuum polarization insertion made in the gluon propagators. 
There will also be additional bremsstrahlung contributions arising from the emission of 
any one of n/ pairs of soft quarks. The calculation of the n/ dependent two-loop ampli- 
tude is technically simplified by writing the amplitude as an integral over a fictitious gluon 
mass. We will now show how this is done for both the virtual loop and bremsstrahlung 
corrections. 

To simplify the virtual correction, consider writing the one-loop virtual contribution, 
in terms of an integral over the virtual gluon momenta, 



where D{k 2 ) is the gluon propagator and w(k;pi) is a weight function depending on the 
gluon momentum, k, and the external momenta, Pi, that can be calculated from ordinary 
Feynman diagrammatic techniques. The one-loop gluon vacuum polarization insertion 
can be made in equation (|2|) by replacing the gluon propagator with D(k 2 )V r (k 2 ), where 
V r (k 2 ), is the dimensionless gluon vacuum polarization renormalized at k 2 = —m 2 . It 
is defined in terms of its unrenormalized counterpart, V(k 2 ), by a subtraction at the 
Euclidean point k 2 = —m 2 , 



Renormalizing in this way corresponds to the Brodsky - Lepage - Mackenzie so-called 
V scheme (in which a s (Q) is normalized by the Coulomb potential between two infinitely 
massive color objects at momentum transfer Q). The n/ dependent part of the two- loop 




(1) 




(2) 



V r {k 2 ) = V(k 2 ) - V{-m 2 t ) 



(3) 



virtual correction is, when expressed in the same form as equation (|3|), 
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^vlt = [ d A ka s (m t )w(k;pi) - . (4) 



k 2 

The vacuum polarization of the gluon can be written as an integral over its imaginary 



part through the dispersion relation with a subtraction at k 2 = —m 2 



t ■ 



1 r°° 

V r (k 2 ) = - / ImP(/i 2 

7T JO 



fi 2 — k 2 /i 2 + m 2 



dfi 2 . (5) 
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By adding and subtracting 1/yU 2 in the weight function, the dispersion relation in 
equation (|5|) becomes, 

v r ik 2 ) = -- r ImP(/i2) + ™t r lmV ^ df f. (6 ) 

7T jo n 2 (k 2 — jj?) ix Jo /i 2 (/i 2 + m 2 ) 

The first term in the dispersion integral (0) when taken together with equation (^) 
has the form of an integral over the gluon mass. The denominator in the first integrand 
in equation (^) is a massive gluon propagator. The rif dependent part of the virtual 
two-loop contribution to the width can be rewritten as an integral over the gluon mass in 
the one-loop contribution, 



5r (2) - - 1 



virt _ , 

7T JO 



r« (a) - m t r (1) (o) 

x virtyh" } o | 2 vvrt\ ) 



fi 2 + m 2 



d^ 2 , (7) 



where T^i rt (/j<) is the virtual contribution of the one-loop width calculated with a gluon 
of mass /i. Although T^ t (/i) is divergent as jj 2 — > 0, the second term in the integrand 
of equation (|?p should be understood as having a small regularizing gluon mass. This 
infrared divergence will be cancelled by the bremsstrahlung terms, below. 

We now show that an analogous relation is true for the bremsstrahlung contribution. 
The rif dependent part of the 0(a s ) bremsstrahlung correction receives contributions from 
two sources. The first corresponds to the emission of two soft tag quarks. The second 
comes from making a vacuum polarization insertion in the final state gluon propagator 
of normal one gluon bremsstrahlung, and corresponds to the renormalization of the gluon 
wave function in the process t — > bWg. The former type of process, quark bremsstrahlung, 
shall be discussed first. 

The four body phase space can be written as the integral over the product of a three 
body phase space and a two body phase space, so that the first bremsstrahlung contribu- 
tion to the width is given by, 
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r 2) f m t M 2 , , . , , . du 2 

Fq-brem = J Q 7^ d W) ^s(^) (8) 

where M 2 is the square of the matrix element, dr 2 (fi) is the two body phase space of the 
two tag quarks with center of mass energy /z, and dr^mt) is the three body phase space 
of the W boson, the b quark, and an on-shell "gluon" of mass \i. 

The phase space integral can be simplified from the form of Equation (H) . Integration 
over the three body phase space will result in a factor of r^g m (/f), the width of the process 
t — > Wbg(fi), where g(n) is a gluon of mass fi. The integral over the two body phase space 
of the light quarks can be written in terms of the imaginary part of the gluon polarization 
operator. The quark bremsstrahlung contribution is, 

r?L m = — I rl(/i)imP(^)^ 2 • (9) 

* TT JO 

(It should be noted in this connection that in the standard definition of the vacuum 
polarization V(k 2 ) as entering the exact propag ator as D(k 2 ) = (k 2 (1 - V{k 2 ))Y\ the 
contribution of a quark pair to the imaginary part IvcvV is negative.) 

The final contribution comes from the gluon wave function renormalization in the 
single gluon bremsstrahlung. The rif dependent piece is calculated by making a vacuum 
polarization insertion in the final gluon propagator appearing in the soft bremsstrahlung 
process t bWg. The correction to the width is r^ em (0)7 :, (0) and can be written using 
the dispersion integral (H) evaluated at k 2 = 0. 

fi^g-brem = ~ [ F^ em ImV(fl 2 ) 2 ( , 2 j_ ^2\^^' ^ (^) 

where T^ em (0) is the tree- level single gluon bremsstrahlung partial width. 

Like the virtual contribution, equation (0), the quark bremsstrahlung contribution, 
equation (|9|), has the form of an integral over the gluon mass. Note that the upper limit 
of the integration in equation (|9|) is m 2 instead of infinity. However, if the decay rate, 
r^g m (/i 2 ), is understood as containing a step function, i.e. r^ em (/i 2 ) = for all fi 2 > m 2 , 
the total rif dependent two-loop contribution becomes a single integral over the gluon 
mass, 

w (2) = -- r l (rgLOO + r&00 - rr^ r(1) (o)) i^ 2 )^ (") 

7r Jo \ (/r + mf) J \l a 

While individual contributions may be infrared divergent, the integral in equation fllTD is 
well behaved both in the infrared and ultraviolet regions. 
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The procedure for calculating the n/ dependent correction is now clear. First, we 
calculate the leading virtual and bremsstrahlung corrections with a fictitious gluon mass, 
jj,. This gluon mass will also serve the purpose of regularizing any infrared divergences that 
might otherwise arise in intermediate results. Next, we perform a weighted integration 
over the gluon mass. The weight function is given in equation (|TT|). 

The virtual massive gluon correction to the width, T^ rt {jjL 2 ), is found by calculating 
diagrams corresponding to the dressing of the vertex, the initial state top quark propa- 
gator, and the final state b quark propagator with a massive gluon. The explicit result of 
this calculation is, 



r 2a Sf 5 x(28 + 10a - 5x 2 ) 

1 o [ — {o + 2x) H = arctan 

3 4 2 Jx{A - x) 



^2- r 



~^ logx — 2(1 + x) log 



x 



-(x + yjx(x - 4)) log -(x - \fx{x - A)) 



(12) 



where x is the dimensionless gluon to top mass ratio, /i 2 /m 2 , and T is the tree level top 
width. 

The bremsstrahlung term is calculated by considering the emission of a massive gluon 
from either of the quarks in t decay. The partial width of this process is, 

2 



<l) (u) 

brem \r"J 

x{20 + 



r o| — [t(25 - 18x - 7x 2 ) - 6(1 + x) I arctan 
3 Ti 4 



-\/3 + Jx(4 - x) 



x 



+ 



2x 



x 



arctan 



2Jx(A - x) 



(1 



x(4 



x(x — 3) 



+ J (10 + x 2 ) logx + ^(1 + x)(loga;) 2 . (13) 



The sum of the expressions ( |T2"D and (|TB|) is finite both in the infrared and ultraviolet 
regions. It is a simple matter to carry out the integration numerically. For a single light 
tag fermion in an Abelian theory, the imaginary part of the gluon vacuum polarization, 
V(fi 2 ) is —a/3. The proper coefficient for QCD can be obtained by replacing lmV(k 2 ) 
with a s b/4:, where b = 11 — |n/ is the first coefficient in the QCD beta function. The 
integral in equation (0) can be evaluated numerically, which gives the decay rate as 



r r 



3tt 



a s (m t ) 



, 3 2 



\ + ^2.54 N 

47T 



(14) 



Within the BLM prescription, the appropriate normalization point for a s is chosen so 
that there are no rif dependent terms of 0(a 2 ). The normalization point can be found 
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by rescaling ot s equation ([TJ|) by, 



b Q 2 
1 + —a s (Q) log^ 



An 



(15) 



and choosing a scale Q such that the rif dependent term in equation fll4|) is cancelled 
by the rif dependent term coming from the rescaling of a s . Numerically, this scale is 
0.281m t . Equation (14) can be written with this new scale, 



1 - 



2tt 2 



3tt 



a^(0.281m t ) 



(16) 



The coupling constant in the above equation is written in terms of Oig , which is derived 
from the potential between two infinitely massive quarks. This coupling constant can be 
related to the more conventional MS scheme by aY(k) = af IS (e~ 5 ^ 6 k) @. The coupling 
constant in equation (|T6D is then equal to af IS (0.122m t ). 

It is interesting to note that this scale is close to the scale of 0.154m 4 which was shown 
in a previous workB to characterize QCD correction to the electroweak p parameter. We 
recall the result that the p parameter is, with a properly normalized one gluon correction, 



5p = 



3 G F m\ 



l-|^(y + l)«f S (0.154m t ) 



(17) 



to leading order in vn? t jvr? w . 

QCD corrections to the top width and the p parameter are physically characterized by 
momenta transfer scales on the order of m t . When these physical quantities are written in 
terms of a mass defined at long distances (like the pole mass used above), corrections are 
introduced corresponding to gluon effects on momenta scales smaller than m t . This has 
the effect of bringing the normalization point of the coupling constant down to a lower 
scale. As discussed by Sirlini, the QCD correction in equation fll7D can be rewritten 



in terms of the running MS mass, m^ s , normalized at the scale m t . The corrections to 



the p parameter, when expressed through this running MS mass, contain only quantities 
defined at short distances, and any long distance nonperturbative effects are suppressed 
by a factor of 0((A QCD /m t ) 4 ). 



While this form may be theoretically more palatable, the quantity mf is not easily 
interpreted in a physical way. The numerical closeness of the QCD corrections to T t and 
Ap suggests that it may be prudent to express Ap in terms of the top quark decay rate. 
Equations ( |T6D and flTTD can be combined in the relation, 



Ap = ^- 

m t TT 
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The order a s correction is numerically very small (the relative correction is 0.0446a s ). 
Writing the p parameter in this way is convenient because the leading QCD corrections 
are tiny, the normalization point of a s is appropriate for momenta transfers on the order 
of m t . 

Both the p parameter and the ratio T t /m t are physical quantities which only contain 
corrections that enter through short distance physics. If these short distance corrections 
are written in terms of quantities defined at longer distances like the pole mass, there is 
some uncertainty of 0{KQCD/f^t) associated with relating a long range parameter to a 
short range one. This uncertainty in minimized by writing short distance parameters in 
terms of each other, like writing the p parameter in terms of T t /m t . 

This work was supported in part by DOE grant DOE-AC-02-83ER40105. BHS is 
supported by the University of Minnesota Doctoral Dissertation Fellowship program. 
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